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Abstract 

We study the spectral density function of a two-dimensional electron liquid 
in a weak magnetic field, the filling factor ^ 1. A hydrodynamic model for 
low-energy excitations of the liquid is developed. It is found that even at ^ 1 
the density of states exhibits a gap at low energies. Its width 2Eq depends on 
the strength of interaction only logarithmically, 2Eq = {huJc/i')^'a^{i'e'^ /efiVF). 
The effects of temperature and disorder on the density of states are discussed. 

PACS numbers:73.40.Hm, 73.40.Gk 
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Since the discovery of the integer and fractional quantum Hall effects [|l| , the properties 
of two-dimensional (2D) electron systems in a strong magnetic field have attracted persistent 
attention. For a long time, the magnetotransport coefficients were the main objects of study. 
Recently, however, experiments allowing a direct investigation of the one-electron density of 
states were performed In Ref. the tunnel current between two 2D layers at small 

filling factors, u, was strongly suppressed at low bias, implying the existence of a gap in the 
spectral density at low energy. 

This fact, which can only be caused by many-electron effects, motivated extensive the- 
oretical studies 1^-0] of tunneling at small filling factors u < 1. Despite rather different 
approaches, all these papers predicted that the width of the gap should be related to the 
single length scale in the lowest Landau level by Ei ~ e^/ (el) where e is the dielectric con- 
stant. The length scale / can be taken to be either the interparticle spacing, the magnetic 
length, or the cyclotron radius since these are approximately equal in this regime. Such a 
value of the gap can be easily understood: Ei equals the Coulomb energy of interaction of 
the extra electron with a "frozen" 2D electron system. If there were no relaxation in the 
electron system at all, the density of states would vanish at energies \huj\ ^ Ei. The re- 
laxation processes smear out the threshold in the density of states which therefore becomes 
sensitive to the details of these processes. However, the suppression is still strong at energies 
\huj\ ^ El if the characteristic relaxation time exceeds h/Ei. 

In this Letter we study the energy dependence of the density of states at large non- 
integer filling factors, ue^ /{ehvp) ^ 1 {vp is the Fermi velocity at zero magnetic field). 
In this regime the interparticle spacing, the cyclotron radius, and the magnetic length are 
distinctly different from each other. We will show that a gap followed by a sharp peak exists 
in the one-electron density of states. Furthermore, in contrast to the case z/ ^ 1 the width 
of the gap depends only logarithmically on the strength of interaction between electrons. At 
u ~ ehvp/e^i which is the lower boundary of the applicability region, our result approaches 
the value Ei ~ /el with the length scale / equal to the cyclotron radius. 

Qualitative Discussion — Our basic approximation is to treat the electrons as a charged 
continuous liquid. Immediately after the electron tunnels into the electron liquid, the system 
acquires extra energy of order e^ri}J'^ je due to the interaction of this electron with the ones 
forming the liquid (here rif. is the 2D electron concentration). The energy of repulsion is 
much greater than the energy of the final state which is determined by an external bias. 
Therefore, in order to reach the final state, the perturbation of the net electron density 
has to spread. The spreading occurs through the virtual emission of collective excitations, 
magnetoplasmons, which have a gap equal to the cyclotron frequency lOc in their spectrum: 



where m is an effective electron mass. Because of the gap, the emission of magnetoplasmons 
occurs within a time of order . Therefore the density perturbation spreads over an area 
of radius — VgUJ~^ where Vg is the group velocity of the magnetoplasmon (|l]). Calculating 
the value of Vg from Eq. (p, we find 




(1) 
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As the charge spreads on the time scale ~ ^^c^, the Lorentz force causes circular currents to 
appear. The Lorentz force due to these circulating currents compensates the excess pressure 
in the region of the density perturbation, blocking the further spreading of the charge and 
making this "vortex" configuration stable. 

The total energy of the vortex consists of the Coulomb energy plus the kinetic energy 
of the circulating currents. The former part is of the order of e^(£ac)~^. The use of Eq. 
yields the Coulomb energy of the vortex ~ mujl/rie. Remarkably, this expression does not 
contain the gas parameter at all. The kinetic energy can be estimated as ~ mv^{nf.al), 
where v is the characteristic velocity of the liquid within the vortex. Using the fact that the 
Lorentz force compensates the force of electric repulsion ~ e'^{eal)^^, we can write down 
the characteristic velocity as f ~ / {muocal). Then, Eq. @) shows that the kinetic energy 
approximately equals the Coulomb energy. The more rigorous treatment below reveals a 
logarithmic divergence of the kinetic energy at small distances. This divergence should 
be cut off at distances of order the Larmour radius Rc = vp/ujc since this is the scale of 
the spatial dispersion of the kinetic coefficients, and thus the length scale on which the 
hydrodynamic approximation fails. The final result for the energy of the vortex Eo is 




This qualitative argument requires the conditions Eq <^ and ^ Rc which are both 
met provided the inequality ve^ / [ehvp) ^ 1 holds. 

The value given by Eq. (|^) sets a new energy scale associated with the tunneling electron. 
If the characteristic time of the further spreading of the vortex (caused, e.g., by dissipation) 
is much greater than h/Eo, a sharp peak appears in the spectral density at the energy Eq. 

Electron Liquid Model — It is clear from the above discussion that the density of states 
at energies E Eq is determined by the long wave length excitations of the electron system. 
We describe these excitations by a hydrodynamic model defined by the velocity v{r,t) and 
electron density n{r,t). The deviation of the concentration 6n{r,t) from its equilibrium 
value ne{r) is related to v{r,t) by the linearized continuity equation 6h + V^, {ngv) = 0. 
The equilibrium electron density rie may depend on the coordinate because of a long range 
disorder potential produced by remote impurities. However, the spatial scale of the disorder 
potential is assumed to be larger than R^, so that the hydrodynamic description still holds. 

It is convenient to take account of the continuity constraint by introducing the field of 
liquid "displacements" u{r,t), 

V {r,t) = u{r,t) , Sn{r,t) = -V {ne{r)u{r,t)) . (4) 

The Hamiltonian of the liquid in a magnetic field is 

H = d^r { — p ^ ' \zxu\ \ + 

J [2mne(r) \ 2 ^ 7 

e2 r (rie (r) u (r)) V^, (rie (rQ u (ri)) | 
t/^^^ 1' 

where z is the unit vector perpendicular to the plane of the liquid. The operator p is 
canonically conjugate to the operator u, [pi{r-i),Uj{r2)] = —ih5ij5{r\ — r^) ■ We used the 
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cylindrical gauge for the vector potential in the kinetic energy term in Eq. the last term 
in the Hamiltonian is the density- density interaction. 

Tunneling Formalism — For the quantitative description of tunneling we have to evaluate 
the spectral density function [ p!0[] , 



A{uj) = 2Re / dte'^^'iifj (t) ifj^ (0) + (0) ip (t)), 



(6) 



where ip are electron creation and annihilation operators, respectively, and (...) denotes 
averaging over the Gibbs distribution. We treat electrons as spinless fermions. 

Now we have to establish the relation between the fermionic operators tp"^ , and the 
hydrodynamic variables u and p. We consider only the local density of states, for which 
ip"^ and ip are taken at r = 0, and adopt the following |TI| approximation for the operator 
creating an electron at the origin: 



^^{t) oc exp[ffl(t)//i] 
with the "shift" operator 11 defined as 

d'^r I r • p{r,t) mujc 



Uit) 
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r- [z X u{r, t)] 



(7) 



(8) 



The operator Eq. (|^) creates an electron density perturbation localized near the origin and 
does not excite any motion of the electron liquid: 

^^{t),6n{r,t)] =-S{r)i,^it), [^j^ {t),v {r,t)] = 0, 

where the velocity operator is v{r) = i[H,u{r)]/h. Because we are considering only the 
local density of states, our calculation is strictly applicable to the case of tunneling from a 
point contact to a 2D plane; however, because the magnetic field restricts the extent of the 
electron — a plane-wave basis is not suitable, for instance — we believe that this calculation is 
relevant to the experimentally interesting case of plane to plane tunneling. For the case of the 
local density of states, we expect the representation Eq. to be an accurate description 
of the long-wavelength collective effects responsible for the gap structure in the density 
of states, while the effects of the single-particle degrees of freedom should not affect this 
structure. 

The representation Eq. (|^) of the creation operator ip"^ in terms of the shift operator 
n(t) allows us to relate the spectral density function A{uj) to the response function of the 
electron liquid. Because the Hamiltonian Eq. (|^) is quadratic, and operator (H) is linear in 
u and p, the calculation may be performed in a standard way |]TU[ and results in 



A{uj) oc Re / dte'^' (e^^*) + e'^*^*)) 



(9) 



Jit) = i[D{t) - D{0)]/h. 



Here we introduced the Green function of the electron liquid D(t) = —i{Tfl(t)fl{0))/h where 
T stands for the time ordering. The most convenient way to proceed is to relate D{uj) to 
the retarded Green function D^{t) = -i9{t){[Il{t),tl{0)])/h by the identity |lOl 
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D{lu) = ReD^{uj) + % coth ImD^(o;). (10) 

Using the equations of motion for the variables u and p, one can check that 11 = 
/ (i^r(5n(r)(e^/e|r|). This would allow one to express the Green functions D and in 
terms of corresponding density-density correlation functions |TD[ . However, it is more conve- 



nient to obtain an equation that determines directly the Green function D^. To derive this 
equation, we notice that the retarded Green function coincides with the linear response of 
the average (n(t)) to an external perturbation of the form 5(t)n(t) added to the Hamilto- 
nian Eq. @. We calculate this linear response with the help of the equations of motion for 
u and p. Because the Hamiltonian Eq. (|^) is quadratic the equation for the Green function 
can be found in a closed form. A simple calculation yields D^(t) = \imr->o D^(t,r), where 
D^(t, r) is the solution of the equation 



em [\ oi ) J \r — r\\ j 

at t > 0, and r) = at t < 0. 

Homogeneous Liquid — For a homogeneous liquid (rie = const), Eq. (^) can be easily 
solved by Fourier transformation. Because Eq ^ hujc-, only the asymptotic behavior at 
t 3> cu"^ is needed. Using this solution in Eqs. (0) and (|10|), we obtain 




m = ( ^ 1 ( <^^YT^k^ (12) 



where Oc is defined by Eq. (||), and 



The logarithmic divergence of the integral in Eq. (jT^) is cut off at wavevectors of order R~ 
where the hydrodynamic approximation fails. Substituting Eq. ([T2|) into Eq. (|^) yields the 
spectral density for ksT -C Eq, 

with energy Eq defined by Eq. (|^). The lower the temperature, the sharper the peak in the 
density of states at Tluj = Eq. The gap at smaller energies is due to the blockade of spreading 
of the electron liquid. As we already discussed, spreading results in the formation of a finite- 
size vortex with energy Eq. The peak at negative energy corresponds to an "antivortex" 
formed after an electron tunnels out from the electron liquid. 

Finite temperature broadens the peak in the density of states by introducing some initial 
inhomogeneities in the electron liquid. For example, a vortex carrying some charge 5q can 
be formed spontaneously by a thermal fluctuation; the probablity of such a configuration is 
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proportional to exp{—ei/kBT) with energy = EoiSq/e)"^. Tunneling of an electron into the 
center of this configuration produces a vortex containing charge Sq + e, and therefore having 
energy huj = Eq{1 + 26q/e). These qualitative considerations lead directly to Eq. (|T3|). 

As we already mentioned, the sharp peak in the density of states at non-zero energy 
occurs because the charge spreading is blocked by the magnetic field applied to an ideal, 
homogeneous electron liquid. There are several ways in which this result should be modified 
in the non-ideal case. First, in a "viscous" liquid, a finite longitudinal dc-conductivity 
a XX allows charge spreading, and the gap is washed out. The microscopic theory [|1^ of 
conductivity for a special case z/ = 1/2 allowed He et al. ^ to obtain the function A{ijj) in 
a wide range of energies. Second, at small filling factors the electron liquid may crystalize 
into a Wigner crystal. It was shown by Johansson and Kinaret that the peak in the 
tunneling density of states in this case has a finite width. This occurs solely due to the 
existence of gapless magnetophonon modes in a crystal that are absent in a liquid. Neither 
of these approaches is suitable for the large filling factors considered here. Third, the charge 



may spread because of inhomogeneities caused, e. g., by smooth disorder; we pursue 
this possibility below. 

Inhomogeneous Liquid — It is well known that density inhomogeneities may give rise to 
gapless excitations of an electron liquid in a magnetic field. These gapless modes redistribute 
charge, and thus cause spreading of the vortex and broadening of the peak in the density of 
states. The origin of such excitations is analogous to that of the edge magnetoplasmon fl^ 



propagating along the periphery of the system. 

In order to estimate the possible effect of smooth disorder, we start with the simplest 
model of a constant gradient of the concentration ne(r) in the region of the liquid into which 
the electron tunnels: 

ne(r) = ne + |Vne|x. (15) 

Here the x axis is in the direction of the gradient, and we assume the potential to be smooth 
so that the concentration changes slowly on the scale set by the vortex size, |Vne|ac -C Ue- 
In this case Eq. ( pTD for zero temperature and Uct ^ 1 yields 

r d'k a,e{u{k)t) e— W*-l 
= I {2^f\k\{l + a^\k\) \hu^{k)\ ' 

where is defined by Eq. ([T^) , and the spectrum uj{k) of the magnetoplasmon excitations 
which are associated with the inhomogeneity Eq. (^5|) is given by 



Kt t <^ Td the exponent in Eq. ( piq) can be expanded, giving J{t) = —iEot/h — 
Eft'^ / {2TTTdh) , similar to Eq. (p!2|). The corresponding result for the spectral density function 
at energies \ \hu\ — Eo\ Et is: 



A(uj) oc \ — ^ exp 



E, 



/ I \ 2" 

TVTdEt I \ — Eq\ 



2h \ E, 



(18) 
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Thus, the spreading of the charge through the low- frequency modes, Eq. ([TTD, leads to a 
broadening of the peak similar in form to that caused by temperature in the homogeneous 
case, Eq. ([M])- As in the case of finite temperature, one can view this broadening as due to 
fluctuations in the initial conditions; however, the fluctuations now have a quantum origin 
in that they are produced by the zero-point motion of the low-frequency modes. 

The region of validity of Eq. ([T8|) corresponds to relatively short times t ^ contributing 
to the integral in Eq. (|^). On this time scale the vortex spreads on a distance Id — 
Td{duj / dky)\\k\^ila^ ~ ttc. According to our assumption about the smoothness of the disorder, 
the corresponding spatial variation of electron density is small, which justifies the use of the 
gradient expansion Eq. (p^). 

We have found that the width of the peak in the local density of states is determined by 
the local concentration gradient in the case of smooth disorder. Turning now to large-area 
tunnel barriers, the differential conductance is related to the ensemble average of the density 
of states, and therefore the width of the peak of the average density of states is determined by 
the typical value of the concentration gradient. This yields a relation between the observable 
width of the peak and parameters of a sample. 

For estimates we assume that the 2D electron liquid is created in a heterostructure, 
where the random potential is produced by a layer of the remote charged impurities with 
concentration N£,, and that all the electrons in the 2D liquid originate from these impurities, 
(ng) = Nj:,. Then, the characteristic spatial scale of the disorder potential is equal to the 
spacer width d and the typical deviation of the concentration from its average value (ne) is 
\/ Nd / d. The density fluctuations are small if N^d"^ ^ 1. We can determine the width of 
the peak SE ~ {Eth/rdY^'^ after estimating the characteristic values of Et and r^: 

6E ae ( 
Eq d \ ehvp 

where as = eU^/e^m. As discussed above, this width is caused by quantum fluctuations 
of the low-frequency modes of the inhomogeneous electron liquid. Note that this width is 
significantly larger than that caused by the direct effect of static disorder on Eq [i.e., by 
ne(r) variation in Eq. (Q)]. The constraint of smooth disorder on the scale set by the vortex 
size is satisfied if ^ {ehvp / e'^){d/ aB^^"^ ■ The opposite case of a "short-range" potential 
Rc ^ d <^ ac will be considered elsewhere. 

In conclusion, we have demonstrated that even relatively weak magnetic field induces a 
gap in the tunneling density of states of an electron liquid. For a homogeneous liquid this 
gap is followed by a sharp peak at a finite energy. Low-frequency modes induced by disorder 
broaden this peak through quantum fluctuations. Eq. (|T9|) shows that for a high-quality 
heterostructure {d/as — 10) the peak should be observable for filling factors i/ ~ 3 — 4. 

Discussions with R. C. Ashoori, J. P. Eisenstein, B. I. Halperin, K. A. Matveev and B. 
I. Shklovskii are gratefully acknowledged. The work at the University of Minnesota was 
supported by NSF Grant DMR-91 17341. 
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